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We study singlet and triplet correlation functions of static quark anti-quark pair defined through 
gauge invariant time-like Wilson loops and Polyakov loop correlators in finite temperature SU(2) 
gauge theory. We use the Liischer-Weisz multilevel algorithm, which allows to calculate these corre- 
lators at very low temperatures. We observe that the naive separation of singlet and triplet states 
in general does not hold non-perturbatively, however, is recovered in the limit of small separation 
and the temperature dependence of the corresponding correlators is indeed very different. 
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I. INTRODUCTION 

It is well established that strongly interacting matter 
undergoes a deconfining transition at some temperature 
which is triggered by large increase in number of degrees 
of freedom (e.g. large increase in the entropy density) as 
well as melting of hadronic degrees of freedom. One of 
the most important feature of the deconfined phase is the 
screening of color charges. It has been argued that color 
screening will lead to quarkonium dissociation above the 
deconfinement temperature which can be used as a sig- 
nature of quark gluon plasma formation in heavy ion col- 
lisions Melting of quarkonium states can be rigor- 
ously studied in terms of spectral functions. Attempts 
to reconstruct spectral functions from Euclidean time 
quarkonium correlators calculated on the lattice have 
been presented in Refs. @, I 1 (for light mesons see 
@) S 13 ) • These studies seemed to indicate that charmo- 
nium states may survive up to unexpectedly high tem- 
peratures 1.6T C — 2.2T C (T c being the transition temper- 
ature). However, it turns out that reconstruction of the 
spectral function is difficult Q . The only statement that 
can be made with confidence is that the Euclidean time 
quarkonium correlation functions do not show significant 
temperature dependence Q. 

On the lattice color screening is usually studied in 
terms of the Polyakov loop correlator related to the free 
energy of static quark anti-quark pair Q . Unlike quarko- 
nium correlators it shows very significant temperature 
dependence across the deconfinement transition. In fact, 
in the deconfined phase the free energy of static quark 
anti-quark pair shows large temperature dependence even 
for very small separations between the static quark and 
anti-quark, much smaller than the inverse temperature 
[lol [ill [l2[ ■ In perturbative picture this can be under- 
stood due to the fact that in the deconfined phase not 
only singlet quark anti-quark (QQ) states contribute to 
the free energy but also colored states with QQ in the 



adjoint (octet for SU(3) and triplet for SU(2)) repre- 
sentation. This observation is also supported by lattice 
calculations of the correlation function of two temporal 
Wilson lines in Coulomb gauge, which in perturbation 
theory corresponds to the so-called singlet free energy 
[H M, EE M, M, GJ, 03 ■ The singlet free energy is 
temperature independent at short distances and coin- 
cides with the zero temperature potential as expected. 
However, at larger distances, e.g. distances of the order 
of typical quarkonium size, it also shows significant tem- 
perature dependence. Thus there seems to be a puzzle: 
heavy quarkonium correlators show almost no tempera- 
ture dependence, while static mesons are largely affected 
by the deconfined medium. 

In the past several years quarkonium properties at fi- 
nite temperature have been studied in potential models 
which use the singlet free energy as an input [13, [U, 
HI HI HI H! HI- Furthermore, within potential model 
calculations it has been shown that the small temper- 
ature dependence of the quarkonium spectral function 
does not necessarily i mply survival of quarkonium states 
at high temperatures [25|, Hl[ . Most quarkonium states 
are melted due to color screening as it was originally sug- 
gested by Matsui and Satz but threshold enhancement 
can compensate for absence of bound states (2f| HH . 

Potential models can be justified in the effective field 
theory framework, the potential non-relativistic QCD 
(pNRQCD), where both scales related to the heavy quark 
mass and bound state size are integrated out [27]. This 
approach can be generalized to finite temperature, where 
also the thermal scales (the temperature and the Debye 
mass) have to be eventually integrated out (28|. In pN- 
RQCD both color singlet and color octet QQ states are 
present as an effective degrees of freedom. The prob- 
lem of defining color singlet and adjoint QQ states on 
the lattice has been considered in Ref. [2{|. It has been 
found that the conventional definition of singlet and ad- 
joint states have problems. Since this question is very 
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relevant for quarkonium physics at finite temperature as 
well as from purely conceptual point of view a more de- 
tailed study is needed. 

In this paper we study static meson correlators in 4 
dimensional SU(2) gauge theory at finite temperature 
and show how the problem observed in Ref. [29J can be 
resolved in the limit of small distances and/or high tem- 
peratures. The rest of the paper is organized as follows. 
In section II we discuss static meson correlators at finite 
temperature and their interpretation in the limit of high 
and low temperatures. Section III contains our numeri- 
cal results along with some important technical details. 
Finally in section IV we give our conclusions. 



II. STATIC MESON CORRELATORS 

Following Ref. [2!| we start our discussion consider- 
ing static meson operators in color singlet and adjoint 
(triplet) states 



0(x, y; t) = i/)(x, t)U(x, y; t)ip(y, t) 



(1) 



O a (x,y-t) = $(x,t)U(x,xo;t)T a U(x ,y;i)iP(y,t),(2) 

where T a is the group generator, ip and ip are creation 
and annihilation operators of static quarks and U (x, y; t) 
are the spatial gauge transporters. In this section we 
consider the general case of SU (N) group, although the 
numerical calculations have been done for N — 2. Next 
we consider correlators of these static meson operators at 
time t — 1/T which after integrating out the static fields 
have the form [29]: 

d(r, T) = ^(TtL\x)U(x, y; 0)L{y)U\x, y, 1/T)),(3) 



G a (r,T) = 



1 



N 2 - 1 



(Trtf(x)Tr L(y)) 



N(N 2 - 1) 



{TiL\x)U(x,y; 0)L(y)U*(x, y,l/T)),(4) 
r=\x- y\. 

Here L(x) is the temporal Wilson line, which on the lat- 
tice is simply L(x) — IXr=o Uo(x, t) with Uo(x,t) be- 
ing the temporal links. The correlators depend on the 
choice of the spatial transporters U(x,y;t). Typically a 
straight line connecting points x and y is used as a path 
in the gauge transporters, i.e. one deals with time-like 
rectangular cyclic Wilson loops, i.e. Wilson loops wrap- 
ping around the time direction. This object has been 
calculated at finite temperature in hard thermal loop 
(HTL) petrurbation theory in context of re-summed per- 
turbative calculations of quarkonium spectral functions 
[3fl l3ll . [HJ ■ This is one of the reasons we are interested 
in non-perturbative evaluation of it. In the special gauge, 
where U(x,y;t) = 1 the above correlators give standard 
definition of the so-called singlet and adjoint free energies 

exp(-F 1 (r,T)/T) = ^(TvL^x)L(y)) , (5) 



exp(-F a {r, T)/T) 
1 



1 



N 2 - 1 



(Trit(x)Tri(y)) 



N(N*-l) {TcL *V> L W)- 



(6) 



The singlet and triplet free energies can be calculated at 
high temperature in leading order HTL approximation 
[33j | resulting in: 



F 1 (r,T) = - 



N 2 



2N 



la s (N 2 - l)a s m D 

CXp( — TYl JjTJ 



27V 



F a {r,T) = +— — exp(-m D r) 



(N 2 - l)a s m D 
2N 



(7) 

(8) 



with mo = yT\J {N/3) being the leading order Debye 
mass. At this order F\ and F a are gauge independent or 
in other words do not depend on the choice of the par- 
allel transporters U(x, y; t). Note that at small distances 
(rmjj <^ 1) the singlet free energy 



fi(r,T)~- 



N 2 - 1 a s 
2N r 



(9) 



is temperature independent and coincides with the zero 
temperature potential, while the adjoint free energy 



F a (r,T) 



1 a s 
2NV 



N 



-a s m D 



(10) 



depends on the temperature. 

The physical free energy of a static quark anti-quark 
pair, i.e. the one related to the work that has to be done 
to separate the static charges by certain distance is given 
by the thermal average of the singlet and adjoint free 
energies 

exp(-F(r,T)/T) = 

1 N 2 - 1 

exp(-F 1 (r, T)/T) + exp(-F a (r, T)/T) 

= ±{TrL{x)TrL(y)) = ±G(r,T). (11) 

This quantity is explicitly gauge independent. In leading 
order HTL approximation the free energy is 



F{r,T) 



(N 2 



AN' 



-^^|exp(~2m D r). (12) 



r 2 T 



The 1/r 2 behavior is due to partial cancellation between 
the singlet and adjoint contribution [34[ and has been 
confirmed by lattice calculations in the intermediate dis- 
tance regime above deconfinement [l2|, HH ■ 

Using the transfer matrix one can show that in the 
confined phase t 29j 



Gi(r,T)=E^i^Me- B " ( ^ T)/T , 
G(r,T) = £r=ie- E " (r ' T)/T > 



(13) 
(14) 
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where E n are the energy levels of static quark and anti- 
quark pair. The coefficients c„(r) depend on the choice 
of U entering the static meson operator O in Eq. (1). 
Since the color averaged correlator G(r, T) corresponds 
to a measurable quantity (at least in principle) it does 
not contain c n . The lowest energy level is the usual static 
quark anti-quark potential, while the higher energy lev- 
els correspond to hybrid potentials [35J, [36j, |37|, |38| . Using 
multi-pole expansion in pNRQCD one can show that at 
short distances the hybrid potential corresponds to the 
adjoint potential up to non-perturbative constants [27| . 
Indeed, lattice calculations of the hybrid potentials indi- 
cate a repulsive short distance part (3f| HI, H3, HI] . If 
ci = 1 the dominant contribution to G a would be the 
first excited state E%, i.e. the lowest hybrid potential 
which at short distances is related to the adjoint poten- 
tial. In this sense G a is related to static mesons with 
QQ in adjoint state. Numerical calculations show, how- 
ever, that Ci(r) =/= 1 and depends on the separation r. 
Thus G a also receives contribution from E\ 29]. The 
lattice data seem to suggest that c\ approaches unity at 
short distances [2!| in accord with expectations based on 
perturbation theory, where c\ — 1 up to 0(g 6 ) correc- 
tions [27]. Therefore at short distances, r < 1/T the 
color singlet and color averaged free energy are related 
F(r, T) = F x {r,T) + T\n(N 2 - 1). This relation is in- 
deed confirmed by lattice calculations [14] . In the next 
section we are going to study the singlet and averaged 
correlators G\ and G in the confined phase and extract 
the coefficients c\. 



III. NUMERICAL RESULTS 

We have calculated correlation functions of static 
mesons Gi(r,T) and G(r,T) both in the confined and 
deconfined phase of SU(2) gauge theory. In our calcu- 
lations we used the Liischer-Weisz algorithm for noise 
reduction [39j . which makes possible to calculate these 
correlators at low values of temperature T not accessi- 
ble by standard methods. Calculations have been done 
for (3 — A/g 2 — 2.5 and 2.7 using lattices with spatial 
extent N s — 24 and N s — 32. Gauge configurations 
have been generated using combination of heat-bath and 
over-relaxation algorithms. Measurements are taken af- 
ter 2 complete updates. A complete update consists of 
one heat-bath and two over-relaxation updates of the en- 
tire lattice and hundred slice updates, in which all links 
inside a slice N% x 2 except of the boundary are updated 
for all N T /2 slices with heat-bath. 

We have studied the color singlet and averaged cor- 
relators given by Eqs. |3]) and (fTTj) . The spatial links 
entering the tr ansp orter U(x,y;0) were smeared using 
APE smearing [40J, which has been applied iteratively. 
The weight of the staple in the APE smeared link was 
0.12. For (3 = 2.5 we use spatial links with 10 steps of 
APE smearing and unsmeared spatial links. For (3 = 2.7 
we used unsmeared spatial links as well as spatial links 



with 10 steps and 20 steps of APE smearing. The lattice 
spacing has been set using the string tension calculated in 
Ref. [H| a?a = 2 .5 = 0.0363(3) anda 2 ^ =2 . 7 = 0.0112(2). 
When quoting the results in terms of reduced tempera- 
ture T/T c we use the value T c /^ = 0.69(2) [Uj. The 
simulation parameters, including the different levels of 
APE smearing used in the present study are summarized 
in Table [J 

For comparison with the study of static meson cor- 
relation functions in Coulomb gauge 15| we have also 
performed calculations on 16 3 x 4 lattice at (3 = 2.3533 
and 2.4215. The two gauge coupling correspond to tem- 
peratures 1.2T C and 1.5T C respectively (see Ref. [l5[ for 
details). Here we used 5, 10 and 20 APE smearings. The 
results from these calculations will be discussed in section 

Imcl 



A. Color averaged correlator in the confined phase 

The color averaged correlator has been calculated in 
the confined phase in the temperature interval 0.32T C — 
0.95T C (3 = 2.5 and 0.49T C - 0.98T C for [3 = 2.7. Thus 
we study this correlator also at very low temperatures, 
where they have not been calculated before. The nu- 
merical results for the color averaged free energy for 
(3 = 2.5 are shown in Figure [1] To eliminate the triv- 
ial temperature dependence due to the color trace nor- 
malization in Figure [1] we show the subtracted free en- 
ergy F'(r,T) = F(r,T) — Tin 4, see the discussion in 
the previous section. In the figure also the zero temper- 
ature potential calculated in Ref. [38[ is shown. The 
color averaged free energy does not show any tempera- 
ture dependence up to temperatures of about 0.76T C . For 
T = 0.76T C the color averaged free energy is below the 
zero temperature potential, indicating a decrease in the 
effective strin g te nsion in accordance with earlier studies 
(see e.g. [13, lig). We see even larger change closer to 
the transition temperature, namely for T > 0.95T C . 

Since the temperature dependence for T < 0.76T C 
is relatively small we attempted to fit the color aver- 
aged correlator with the 1-exponential form G(r, T) = 
c"(r) exp(—Ei(r)/T). The ground state energy E\(r) ex- 
tracted from this fit agrees well with the zero tempera- 
ture potential calculated in Ref. [38[, while the coeffi- 
cients c\(r) are close to one as expected (see the discus- 
sion in the previous section). The fit details are shown 
in the Appendix. Although the deviations of c"(r) from 
unity are small, they appear to be statistically significant. 
These deviations increase with increasing r. Therefore 
they are likely to be due to the contribution from ex- 
cited states, as the gap between the ground state and ex- 
cited states (hybrid potential) gets smaller with increas- 
ing separation. Therefore we attempted to fit the color 
averaged correlator with 2-exponcntial form G(r, T) = 
exp{~E 1 (r)/T) + exp{~E 2 (r)/T) using T < 0.76T C . The 
results of the fit are given in the Appendix. We find that 
the value of E\ extracted from the 2-exponential fits are 
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FIG. 1: The color averaged free energy defined by Eq. 
below deconfinement temperature at (3 — 2.5 calculated on 
32 3 x N T lattices. Also shown is the T — potential. 



consistent with the ones obtained from 1-exponential fit 
as well as with the value of the zero temperature poten- 
tial. For E2 (r) we find values which are somewhat below 
the first hybrid potential. This is presumably due to the 
fact that there is small contribution from the higher ex- 
cited states. 



B. Color singlet correlators in the confined phase 

As mentioned above the color singlet correlators have 
been calculated using different levels of APE smearing 
in the spatial gauge connection. It is well known that 
smearing increases the overlap with the ground state by 
removing the short distance fluctuation in the spatial 
links (see e.g. [42|). For this reason smearing also re- 
duces the breaking of the rotational invariance to the 
level expected in the free theory. We have found that 
when no smearing is used the color singlet free energy, 
—ThxGi(r,T) shows a small but visible temperature de- 
pendence. In particular Fi(r, T) is larger than the T = 
potential for intermediate distances 0.5 < r^fa < 2. A 
similar effect has been observed in calculation in 3 di- 
mensional SU(2) gauge theory [2^, [43[ as well in 4 di- 
mensional SU(2) and SU(3) gauge theory calculations 
in Coulomb gauge [HI, |44[ • The temperature dependence 
of the singlet free energy is significantly reduced when 
APE smearing is applied. In Figure [2] we show the color 
singlet free energy for (3 = 2.5 and 10 APE smearings. 
As one can see from the figure the color singlet free en- 
ergy shows significantly smaller temperature dependence 
as we get closer to the deconfinement temperature. In 
particular, only for T = 0.95T C we see significant temper- 
ature dependence, which, however, is much smaller than 
for color averaged free energy. To understand the tem- 
perature dependence of the color singlet correlator we use 
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TABLE I: Simulation parameters, # meas refers to the actual 
number of measurements. 
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FIG. 2: The color singlet free energy below deconfinement 
temperature at (3 = 2.5 calculated on 32 3 x N T lattices. Also 
shown is the T = potential. The inset shows the color sin- 
glet free energy from which the contribution from the matrix 
element Tlnci has been subtracted. 
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1-exponential fit G x {r,T) = ci(r)exp(-Ex(r)/T). In all 
cases considered, the value of £a(r) extracted from fits 
are in good agreement with the calculation of the zero 
temperature potential in Ref . [38] . The value of the pref- 
actor Ci (r) is shown in Figure [3] When no APE smear- 
ing is used the value of ci(r) strongly depends on the 
separation r. At small distances it shows a tendency of 
approaching unity as one would expect in perturbation 
theory. However, c\(r) decreases with increasing distance 
r. At large distance its value is around 0.3 — 0.5 (see Ap- 
pendix for details). Similar results for Ci(r) have been 
obtained in study of SU(2) gauge theory in 3 dimensions 
[29l ] . When APE smearing is applied the r-dependence 
of the amplitude ci (r) is largely reduced and its value is 
close to unity both for (3 = 2.5 and (3 = 2.7. For (3 = 2.7 
we also see that increasing the number of smearing steps 
from 10 to 20 reduces the deviation of c± (r) from unity. 

As discussed in section II perturbation theory predicts 
that the deviations of Ci(r) from unity is of order a^. 
Therefore it can be made arbitrarily small by going to 
sufficiently small distances but even for distances accessi- 
ble in this study these deviations are expected to be small 
based on perturbation theory. It is known, however, 
that lattice perturbation theory converges very poorly. 
The main reason for this has been identified with the 
short distance fluctuations of link variables, which makes 
their mean value very different from unity [451 ]. Smear- 
ing removes these short distance fluctuations and this is 
the reason why ci (r) is much closer to unity when APE 
smearing is applied. 

From the above discussions it is clear that almost the 
entire temperature dependence of the singlet free energy 
at distances 0.5 < r^fa < 2 is due to the deviation of ci 
from unity and can be largely reduced by applying APE 
smearing to the links in the spatial gauge connections. 
To further demonstrate this point in the inset of Figure 
[5] we show the results for F\ (r, T) + T In C\ (r) . Clearly 
no temperature dependence can be seen in this quantity 
up to 0.95T C , where we see temperature dependence at 
distances r^fa > 1.5 corresponding the expected drop of 
the effective string tension. 



Color singlet free energy in the deconflned 
phase 
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FIG. 3: The pre-exponential factor of the color singlet corre- 
lators as function of distance r for /3 = 2.5 (top) and j3 = 2.7 
(bottom). Shown are results for unsmeared spatial links and 
10 and 20 steps of APE smearing. 



Its value is therefore independent of the definition of the 
singlet correlator G\ (r, T) and is related to the renormal- 
ized Polyakov loop L ren (T) = exp(-F oc (T)/(2T)) Q. 

At leading order Fi (r, T) — Foo (T) is of Yukawa form 
(c.f. Eq. 0). Therefore we find it useful to show our 
numerical results in terms of the screening function 



In this subsection we discuss the properties of the color 
singlet free energy, F\{r, T) — —T\nGi(r,T) above the 
deconfincment temperature. It turns out that the singlet 
free energy calculated from cyclic Wilson loops shares the 
same qualitative features as the sin glet free energy calcu- 
lated in Coulomb gauge 0, EE, E 113, QJ] ■ At short dis- 
tances it is temperature independent and coincides with 
the zero temperature potential. At large distances it 
approaches a constant Foo(T) which monotonically de- 
creases with the temperature. The constant Foo(T) is 
the free energy of two isolated static quarks, or equiva- 
lently of a quark anti-quark pair at infinite separation. 



S(r,T) = r-(Fi(r,T)--Foo(T)) (15) 

In Figure 3] we show the results for the screening function 
at different temperatures. At short distances (rT < 0.5) 
the singlet free energy does not depend on the smearing 
level. Furthermore, it is very close to the free energy 
calculated in Coulomb gauge. We expect that at large 
distances the screening function S(r, T) will show an ex- 
ponential decay determined by a temperature dependent 
screening mass mi(T), which is equal to the leading or- 
der Debye mass up to the non-perturbative g 2 correc- 
tions: mi = mn + 0{g 2 ) [i^ . From Fig. 0] we can see 
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FIG. 4: The screening function S(r,T) = r(Fi(r,T)-Foo(T)) 
at different temperatures calculated for /3 = 2.7 as function 
of the distance in units of \[5 (top) and as function of rT 
and different values of f3 (bottom). In the lower panel we 
also shown the results from calculations in Coulomb gauge 
[TEj as open symbols. In the inset the screening masses mi 
extracted from singlet free energies are shown together with 
the results obtained in Coulomb gauge 15] (open symbols). 
The line shows the leading order results for the Debye mass. 
In the upper panel the filled symbols correspond to 20 APE 
smearings steps, while the open ones to 10 APE smearings 
steps. 



that indeed S(r, T) behaves exponentially with screening 
mass proportional to the temperature. We note, however, 
that there is some dependence on the smearing level at 
larger distances. This disappears at high temperatures 
and with increasing the smearing level. In particular, 
for (3 < 2.5 it turns out that there is no dependence 
on smearing level for 5 or more smearing steps. For 
(3 = 2.7 we need 10-20 steps depending on the tempera- 
ture to achieve stable results. Fitting the large distance 
behavior of the screening function by an exponential form 
exp(— mi(T)r) we determine the screening mass m\(T). 
Typically we considered distances r > 1/T as well as the 



maximal number of smearings steps (10 for (3 < 2.5 and 
20 for (3 — 2.7 ) in our fits. In the inset of Fig. d]we show 
the color singlet screening masses extracted from the fits 
and compare them to the results obtained in Coulomb 
gauge in Ref. 15]. In the inset we also show the lead- 
ing order Debye mass calculated using 2-loop gauge cou- 
pling g(/j, = 2ttT) in MS'-scheme. We used the value 
T c /A-g-g — 1.09 [43] to calculate the coupling g. As we 
see from the figure the screening masses are smaller than 
those calculated in Coulomb gauge and agree well with 
the leading order perturbative prediction. 

Leading order perturbation theory also predicts that 
F}(r,T)/F(r,T) = 6 for r > 1/T (c.f. Eqs. [THE]). We 
find that at the highest temperature, T = 3.42T C this 
relation indeed holds within statistical errors. 



D. Color triplet free energy 

We have calculated the color triplet correlator defined 
by Eq. (0J for different temperatures below and above 
the transition temperature. Below the deconfincmcnt 
temperature we observe a moderate T-dependence of the 
triplet correlator. We also find that the corresponding 
free energy — T In G^(r, T) is smaller than the first hy- 
brid potential calculated in Ref. [38| , but larger than the 
triplet free energy in Coulomb gauge [l5j . 

Let us assume that only two states contribute to the 
Eqs. CTJ]) and (H]). Then from Eq. g]) it follows that 

4 F 3 (r, T) = E 2 (r) - T ln(l - c 2 (r) + ~ (1 - c x (r))e AE ^/ T ), 

(16) 

with AE(r) — E 2 (r) — Ei(r). We have seen in sec- 
tion III.B that the temperature dependence of the sin- 
glet free energy is quite small. In any case it is consid- 
erably smaller than the temperature dependence of the 
averaged free energy. Therefore the contribution of the 
excited states to G\ (r, T) is quite small and it is reason- 
able to assume that c 2 {r) <C 1. We also expect that at 
small distances, c 2 (r) ~ (tAqcd) 4 [HI- Thus, the tem- 
perature dependence of Fs(r, T) and its deviation from 
the hybrid state E 2 (r) is due to small deviation of C\ (r) 
from unity. At low temperatures, when AE ^> T these 
small deviations are amplified by the exponential fac- 
tor. To verify this, we have subtracted the correction 
Tln(l + -|(1 — ci)e AB / T ) from the triplet free energy as- 
suming that Ei (r) is given by the ground state potential 
and E 2 (r) is given by the first hybrid potential as calcu- 
lated in Ref. [38| . The numerical results are summarized 
in Fig. [5] which shows that after this correction is ac- 
counted for in the confined phase the triplet free energy 
at low temperatures agrees reasonably well with the first 
hybrid potential. As temperature increases more excited 
states contribute. In particular, at 0.76T C the value of 
the triplet free energy can be accounted for by including 
the next hybrid state [HI]. However, at 0.95T C we see 
large temperature effects, which cannot be explained by 
including the contribution from only few excited states. 



7 




3.5 - 



0.5 1 1.5 2 2.5 3 

r<r 1/2 

FIG. 5: The triplet free energy at different temperatures cal- 
culated at P = 2.5. The filled symbols correspond to cal- 
culations in Coulomb gauge. Also shown is the first hybrid 
potential calculated in Ref. [3Sl |. 

In Fig. [5] we also show the triplet free energy above 
the deconfinement temperature compared to the calcu- 
lations in Coulomb gauge [l5j]. It turns out to be much 
smaller than in the confined phase and agrees well with 
Coulomb gauge results. This means that the small de- 
viation of the overlap factor Ci(r) from unity are unim- 
portant in this case. The triplet free energy monotoni- 
cally decreases with increasing temperature as expected 
in HTL perturbation theory (c.f. Eq. ([8])). In the limit 
of high temperatures and short distances r < l/T we 
have E 2 (r) = a s /(4r), AE(r) — a s /r, c 2 (r) ~ and 
ci(r) = 1 + 0(al). Therefore we can expand the loga- 
rithm in Eq. (fl6|) to get 

F 3 (r,T) = +]^ + 0(a 3 s T) + 0(a s m D ). (17) 
4 r 

Thus the correction due to ci(r) ^ are much smaller 
than the expected leading order thermal effects in the 
triplet free energy. 

While the temperature dependence of F 3 (r, T) and 
Fi(r,T) is different even at short distances the r- 
dependence is expected to be similar. In leading order 
perturbation theory we have 

dFi(r,T) 

dr _ _o 
dF 3 (r,T) 

dr 

We have calculated the above ratio numerically from the 
lattice data and have found that for T > 1.71T C the de- 
viation from the expected value (—3) are smaller than 
10%. 

IV. CONCLUSIONS 

In this paper we have studied singlet and triplet static 
quark anti-quark correlators in finite temperature SU (2) 



theory expressed in terms of Polyakov loop correlators 
and cyclic Wilson loops (cf. Eqs ©, ©)■ The study of 
the latter is interesting as it has been used in resummed 
perturbative calculations of quarkonium spectral func- 
tions [3(1 [Hj]. I n leading order and probably next-to- 
leading order of perturbation theory the static correla- 
tors defined by Eq. §5§ and (d]) give the energies of the 
singlet and triplet states respectively, however, this sepa- 
ration does not hold in general case. Due to interactions 
with ultrasoft fields there will be a mixing of singlet and 
triplet states which is proportional to a\(l/r) and (f/i) 4 , 
with n being the ultra-soft scale [27|. In our case the 
ultra-soft scale can be the binding energy, a s /r, Aqcd 
or g 2 T. Therefore it is expected that mixing is quite 
small at sufficiently small distances. We determined the 
mixing between singlet and triplet states in terms of the 
overlap factor c\ (r) . If the overlap factor is unity there 
is no mixing. Our lattice calculations show that c\(r) 
indeed approaches one at small distances. Using itera- 
tive APE smearing the deviation of c\ (r) from unity can 
be largely reduced. Therefore the contribution of singlet 
state to G3 (r, T) appears to be small at temperatures 
close to deconfinement temperature. This contribution 
is also controlled by the non-perturbative gap between 
the singlet and triplet states, i.e. the gap between the 
static potential and the first hybrid potential. 

Our analysis shows that at short distances rT < 1 
the singlet correlator is almost temperature independent, 
while the triplet correlator is largely affected by the de- 
confinement. The temperature dependence of the triplet 
correlators indicate the melting of the non-perturbative 
gap between the singlet and the triplet states above de- 
confinement, which turns out to be consistent with per- 
turbative expectations. Because of the disappearance of 
the non-perturbative gap the small deviations of ci(r) 
from unity play no role in the deconfined phase. In 
particular, the difference between the singlet and triplet 
correlators defined through cyclic Wilson loops and us- 
ing Coulomb gauge turns out to be quite small. This 
finding is important for application of thermal pNRQCD 
discussed in Ref. [2c| to realistic quarkonia and temper- 
atures not very far from the deconfinement temperature. 

There are obvious extensions of the work carried out in 
this paper. One could do analogous calculations in SU (3) 
gauge theory, however this hardly will give qualitatively 
different results. More interesting would be to extend 
the calculations to QCD with dynamical fermions. First 
steps in this direction are being made using improved 
staggered (p4) action [4^| . It will be interesting to calcu- 
late the correlators G\ (r, T) and G3 (r, T) in perturbation 
theory beyond HTL approximation. One practical prob- 
lem in doing these calculations is the fact that the scales 
T and mo are not separated. This problem, however, 
could be possibly solved by using screened perturbation 
theory (50. l5lll52T| . 
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N s = 24 


Ns = 32 


r 


c i 




x 2 
dof 


C X 


Ex 


x 2 
dof 


1 


1.00077(28) 


0.33524(2) 


0.50 


1.00046(28) 


0.33522(1) 


1.48 


2 


1.0049(13) 


0.48423(8) 


0.49 


1.0039(12) 


0.48415(6) 


0.96 


3 


1.0134(45) 


0.56568(24) 


0.64 


1.0115(32) 


0.56550(16) 


0.93 


4 


1.027(10) 


0.62400(61) 


0.73 


1.0237(67) 


0.62365(35) 


0.75 


5 


1.049(22) 


0.6732(14) 


1.23 


1.038(14) 


0.67232(78) 


0.66 


6 


1.067(43) 


0.7171(27) 


1.43 


1.051(29) 


0.7161(17) 


1.10 


7 


1.087(91) 


0.7582(59) 


1.80 


1.023(60) 


0.7542(40) 


1.53 


8 


1.05(16) 


0.793(11) 


2.06 


0.93(11) 


0.7861(83) 


1.69 



TABLE II: Single exponential fit of G(r, T) at /3 = 2.5 using 
N T = 12 - 24 for two spatial sizes N s = 24 and N s = 32. 



r 


c l 


Ej 


x 1 
dof 


1 


1.00065(44) 


0.28580(2) 


0.64 


2 


1.0084(28) 


0.39584(12) 


1.21 


3 


1.0325(90) 


0.44727(36) 


0.90 


4 


1.078(21) 


0.47939(80) 


0.54 


5 


1.143(39) 


0.5038(14) 


0.41 



TABLE III: Single exponential fit of G(r, T) at (3 ■■ 
N t = 20 - 28. 



2.7 using 
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Appendix 

In this appendix we give some details of the fits of 
the overlap factors and energy levels for different lattice 
volumes. 

Let us first discuss the fits of the color averaged cor- 
relator G(r,T). We have fitted our lattice data at 
low temperatures with one exponential form G(r, T) = 
cf(r) exp(— Ei(r)/T). For f3 — 2.5 we used lattices with 
temporal extent N T = 12, 14, 16, 18, 20, 22 and 24. For 
/3 = 2.7 we used N T = 20, 22, 24 and 28. The values of 
fit parameters c° (r) and E\ (r) are shown in Table fll] for 
(3 = 2.5 and Table ITTT1 for (3 = 2.7. As one can see from 
the Tables at f3 = 2.5 the single exponential fit works 
reasonably well up to distance r = 7 for N s = 32 and 
distance r = 6 for N s = 24. At (3 = 2.7 the fits work up 
to distance r = 5 only. We also performed double expo- 
nential fits G{r,T) = exp(-£i(r)/T) + exp(- E 2 (r) /T) 
to check the reliability of the extraction of E\ . The re- 
sults are shown in Tables IIVI and [V] for (3 = 2.5 and 
(3 = 2.7 respectively. Overall the values of E\ extracted 
from 2-exponential fits are in reasonably good agreement 
with the corresponding values extracted from singlet ex- 
ponential fit. 



We fitted the singlet correlators calculated for different 
number of APE smearing steps with single exponential 
(r.T) = ci(r) exp(-£i <r)/T). The results for 





N 


= 24 




N 


— 32 




T 






x 2 
dof 


-C'l 


p„ 


x 1 
dof 


i 


0.335204(5) 


1.011(41) 


0.89 


0.335199(3) 


1.048(50) 


0.84 


2 


0.483974(18) 


0.989(20) 


0.56 


0.483948(9) 


0.998(15) 


0.47 


3 


0.565022(42) 


0.983(17) 


0.44 


0.564946(29) 


0.982(14) 


0.60 


4 


0.62269(12) 


0.984(21) 


0.55 


0.622522(88) 


0.986(19) 


0.80 


5 


0.67088(42) 


0.987(36) 


1.08 


0.67046(19) 


1.003(24) 


0.65 


6 


0.7140(10) 


1.004(57) 


1.35 


0.71357(48) 


1.023(40) 


0.86 


7 


0.7545(24) 


1.016(86) 


1.69 


0.7533(12) 


1.07(10) 


1.43 



TABLE IV: Double exponential fit of G(r, T) at /3 = 2.5 using 
N T = 12 - 24 for two spatial sizes N s = 24 and N s = 32. 



r 


Ex 


E 2 


x 1 

Hot' 


1 


0.285780(1) 


0.700(11) 


0.14 


2 


0.395551(9) 


0.6819(60) 


0.15 


3 


0.446179(27) 


0.6657(45) 


0.11 


4 


0.476943(53) 


0.6532(33) 


0.07 


5 


0.49965(12) 


0.6445(37) 


0.07 


6 


0.51857(29) 


0.6362(48) 


0.11 


7 


0.53545(67) 


0.6297(63) 


0.16 


8 


0.5514(14) 


0.6240(79) 


0.18 



TABLE V: Double exponential fit of G(r,T) at (3 - 
N T = 20 — 28. 
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the fit parameters c\{r) and Ei(r) at f3 = 2.5 are shown 
in Tables IVT1 and IVlTl for zero and 10 APE smearing steps 
respectively. In Tables IVIIH IIXI and |X] we show the fit 
parameters for j3 = 2.7 and zero, 10 and 20 APE smearing 
steps. 





N s = 24 


N s = 32 


r 


Cl 


Ex 


x J 
dof 


Cl 


Ex 


x 2 
dof 


1 


0.89076(34) 


0.33523(2) 


0.36 


0.89043(28) 


0.33521(2) 


1.31 


2 


0.68239(89) 


0.48406(8) 


0.77 


0.68191(79) 


0.48401(6) 


0.70 


3 


0.5113(16) 


0.56549(20) 


0.20 


0.5095(14) 


0.56519(14) 


1.11 


4 


0.4020(33) 


0.62353(53) 


0.19 


0.4035(24) 


0.62360(33) 


0.58 


5 


0.3439(63) 


0.6731(12) 


0.82 


0.3367(46) 


0.67123(86) 


0.34 


6 


0.313(14) 


0.7174(32) 


2.24 


0.312(10) 


0.7173(21) 


1.32 


7 


0.321(32) 


0.7642(74) 


1.93 


0.261(20) 


0.7496(54) 


1.10 


8 


0.311(46) 


0.803(11) 


0.70 


0.240(38) 


0.786(11) 


3.19 


9 


0.35(13) 


0.845(29) 


1.43 


0.228(79) 


0.822(24) 


1.69 



TABLE VI: Single exponential fit of Gx(r,T) at different lat- 
tice separations r at f3 = 2.5 with no APE smearing using 
N T = 12 — 24. 





N a = 24 


N a = 32 


r 


Cl 


Si 


x J 
dot' 


Cl 


Si 


x 2 
dot' 


1.0000 


0.99457(96) 


0.33528(6) 


0.71 


0.99260(50) 


0.33516(3) 


0.53 


1.4142 


0.9926(19) 


0.42930(11) 


0.81 


0.9904(11) 


0.42916(6) 


0.32 


1.7320 


0.9886(26) 


0.47322(16) 


0.38 


0.9864(16) 


0.47310(9) 


0.33 


2.0000 


0.9919(30) 


0.48414(17) 


0.70 


0.9882(18) 


0.48389(10) 


0.19 


2.8284 


0.9810(69) 


0.55824(42) 


0.53 


0.9721(40) 


0.55758(24) 


0.32 


3.0000 


0.9835(63) 


0.56554(36) 


0.50 


0.9751(46) 


0.56494(27) 


0.28 


3.4641 


0.954(10) 


0.59681(68) 


0.72 


0.9487(60) 


0.59643(37) 


0.47 


4.0000 


0.972(13) 


0.62406(82) 


0.33 


0.9503(98) 


0.62243(60) 


0.19 


4.2426 


0.950(16) 


0.6376(11) 


0.33 


0.938(11) 


0.63656(70) 


0.32 


5.0000 


0.953(26) 


0.6732(17) 


0.83 


0.923(18) 


0.6706(12) 


0.40 


5.1961 


0.891(28) 


0.6803(21) 


1.47 


0.870(16) 


0.6787(11) 


0.72 


5.6568 


0.897(40) 


0.7015(32) 


0.94 


0.880(25) 


0.7000(18) 


0.12 


6.0000 


0.933(44) 


0.7185(31) 


0.48 


0.877(35) 


0.7129(26) 


0.54 


6.9282 


0.765(66) 


0.7473(61) 


1.34 


0.807(41) 


0.7508(33) 


0.88 


7.0000 


0.837(68) 


0.7550(61) 


0.62 


0.874(64) 


0.7563(49) 


0.38 


7.0711 


0.902(87) 


0.7643(71) 


0.99 


0.808(52) 


0.7568(44) 


0.29 


8.0000 


0.71(11) 


0.786(11) 


1.28 


0.86(10) 


0.7963(80) 


1.78 



TABLE VII: Single exponential fit of Gi(r,T) at j3 = 2.5 
calculated with 10 steps of APE smearing using N T — 12 — 24. 



r 


Cl 


Si 


x 2 
dot' 


1 


0.89048(47) 


0.28578(2) 


2.72 


2 


0.6898(18) 


0.39557(11) 


2.30 


3 


0.5241(39) 


0.44628(31) 


2.29 


4 


0.4243(73) 


0.47792(72) 


1.66 


5 


0.371(12) 


0.5019(14) 


1.44 


6 


0.362(21) 


0.5238(25) 


0.83 


7 


0.385(33) 


0.5449(36) 


0.36 


8 


0.396(46) 


0.5612(48) 


1.04 


9 


0.497(76) 


0.5838(64) 


0.74 


10 


0.55(11) 


0.5991(89) 


0.59 



TABLE VIII: Single exponential fit of Gi(r,T) at /3 = 2.7 
with no APE smearing using N T — 20 — 28. 



r 


Cl 


Ei 


2 

X 

dot' 


1.0000 


0.9962(12) 


0.28582(5) 


0.02 


1.4142 


0.9952(25) 


0.35819(10) 


1.00 


1.7320 


0.9927(33) 


0.38969(14) 


0.59 


2.0000 


0.9923(42) 


0.39585(18) 


0.18 


2.8284 


0.9798(97) 


0.44311(42) 


1.74 


3.0000 


0.978(11) 


0.44681(45) 


0.59 


3.4641 


0.971(14) 


0.46541(62) 


1.19 


4.0000 


0.968(20) 


0.47854(88) 


0.41 


4.2426 


0.949(24) 


0.4852(11) 


1.08 


5.0000 


0.944(34) 


0.5018(15) 


0.22 


5.1961 


0.949(37) 


0.5073(17) 


0.89 


5.6568 


0.921(48) 


0.5151(22) 


1.02 


6.0000 


0.920(52) 


0.5211(24) 


0.04 


6.9282 


0.939(76) 


0.5398(34) 


0.78 


7.0000 


0.894(78) 


0.5381(37) 


0.20 


7.0711 


0.894(72) 


0.5399(34) 


0.60 


8.0000 


0.849(99) 


0.5522(50) 


0.02 


8.4853 


0.86(11) 


0.5614(57) 


0.30 


8.6603 


0.89(14) 


0.5661(64) 


0.37 


9.0000 


0.78(13) 


0.5640(73) 


0.24 


9.8995 


0.80(18) 


0.5793(94) 


0.07 


10.0000 


0.75(17) 


0.5767(98) 


0.19 


10.3923 


0.74(19) 


0.584(11) 


0.25 



TABLE IX: Single exponential fit of Gi(r,T) at (3 = 2.7 cal- 
culated with 10 steps of APE smearing using N T = 20 — 28. 



r 


Cl 


Ei 


X 2 /dof 


1.0000 


0.9965(15) 


0.28582(6) 


1.21 


1.4142 


0.9965(26) 


0.35809(11) 


1.10 


1.7320 


0.9957(45) 


0.38949(19) 


2.02 


2.0000 


0.9982(53) 


0.39573(22) 


1.29 


2.8284 


0.993(12) 


0.44262(51) 


1.05 


3.0000 


0.993(13) 


0.44639(55) 


1.38 


3.4641 


0.985(21) 


0.46418(88) 


3.86 


4.0000 


0.985(24) 


0.4771(10) 


2.50 


4.2426 


0.979(31) 


0.4839(13) 


0.65 


5.0000 


0.979(38) 


0.4999(16) 


2.64 


5.1962 


0.984(48) 


0.5049(20) 


3.10 


5.6569 


0.992(58) 


0.5139(25) 


0.86 


6.0000 


0.982(63) 


0.5193(27) 


2.34 


6.9282 


1.04(10) 


0.5378(40) 


2.29 


7.0000 


0.983(86) 


0.5360(37) 


2.10 


7.0711 


1.03(11) 


0.5395(44) 


0.86 


8.0000 


1.00(12) 


0.5525(52) 


2.15 


8.4853 


1.11(18) 


0.5640(70) 


1.08 


9.0000 


0.94(15) 


0.5642(69) 


1.67 


10.0000 


0.92(22) 


0.577(10) 


1.82 



TABLE X: Single exponential fit of Gi(r,T) at (3 = 2.7 cal- 
culated for 20 APE smearing using N T = 20 - 28. 



